A low density finite temperature apparent "insulating" phase in 2D systems 



m 
o 
o 

(N 

in 



D 

■4— » 
•4-* 



c 

O 

o 



(N 



(N 
O 

o 

-I— > 



T3 

c 

o 
o 



X 



S. Das Sarma and E. H. Hwang 
Condensed Matter Theory Center, Department of Physics, University of Maryland, College Park, Maryland 20742-4111 

(February 2, 2008) 

We propose that the observed low density "insulating" phase of a 2D semiconductor system, with 
the carrier density being just below (n < n c ) the so-called critical density where the derivative of 
resistivity changes sign at low temperatures (i.e. resistivity p(T) increases with increasing T for 
n > n c whereas it decreases with increasing T for n < n c ), is in fact a "high-temperature" crossover 
version of the same effective metallic phase seen at higher densities (n > n c ). This low density 
(n < ric) finite temperature crossover 2D effective insulating phase is characterized by p(T) with 
power law temperature dependence in contrast to the truly insulating state (occurring at still lower 
densities) whose resistivity increases exponentially with decreasing temperature. 
PACS Number : 71 .30.+h; 73.40. Qv 



The so-called 2D metallic phase, first identified [1] in 
n-Si MOSFETs by Kravchenko and collaborators, is still 
a subject of considerable interest and controversy [2]. 
Although the effectively metallic character (albeit at fi- 
nite temperatures) of this phase is well-established ex- 
perimentally in many 2D systems, it is still not known 
whether a true T — 2D metal exists. The nature of the 
associated 2D metal-insulator-transition (2D MIT) as a 
function of carrier density (n) is also a mystery (par- 
ticularly since a noninteracting disordered 2D system at 
T = is known on firm theoretical grounds not to un- 
dergo any MIT as all states in 2D are localized in the pres- 
ence of finite disorder [3]). There is increasing evidence 
that the putative 2D "metallic" phase, while being a per- 
fectly good "effective" metal down to T — 50 — lOOmK 
where experiments have so far been carried out, is es- 
sentially a finite temperature effective metal with screen- 
ing [4] and related electron-electron interaction effects [5] 
playing important quantitative roles in determining the 
temperature dependent resistivity (and producing the ef- 
fective metallic behavior). It should also be emphasized 
that no reasonable model or theory for a T = strongly- 
correlated 2D metal (stabilized presumably by electron- 
electron interaction in the presence of disorder) has been 
put forward in the literature [2] . (Such a true 2D metal, 
if it exists at T = 0, cannot be a Fermi liquid since the 
corresponding noninteracting 2D system is an insulator 
at T = [3].) The observed 2D metallic phase shows 
strong temperature dependence (dp/dT > 0) which is 
nominally unexpected in a metal at low temperatures 
where phonons freeze out in the low temperature Bloch- 
Griineisen regime. We have argued, on the basis of con- 
crete calculations, that the strong temperature depen- 
dence in the resistivity in the effective metallic phase 
is caused by a temperature dependent effective disor- 
der arising from the strong temperature dependence of 
low density 2D screening [4]. This view point of the 
metallic temperature dependence has received further 
support from recent theoretical calculations of interac- 
tion corrections to 2D resistivity generalizing the screen- 
ing theory [5]. Unlike 3D metals, where qTp/^kp ~ 1 
and T/Tp <C 1 (where qrF, kp, Tp are the Thomas- 



Fermi screening wave vector, the Fermi wave vector, and 
the Fermi temperature, respectively), the semiconductor- 
based 2D systems at low carrier densities could typically 
have qTF/2kp 3> 1 and T/Tp ~ 1, engendering strong 
screening-induced temperature dependence in the resis- 
tivity arising from the strong temperature dependence in 
the screened impurity disorder scattering in contrast to 
3D metals where any temperature dependence of disor- 
der scattering is exponentially suppressed at low temper- 
atures. 

In this paper we extend our screening calculations to 
very low densities (below n c , the apparent critical density 
for the MIT) , where the experimental resistivity increases 
with decreasing temperature (dp/dT < 0) indicating the 
existence of a nominally insulating phase. Quite sur- 
prisingly we predict the existence of an effective "insu- 
lating" phase defined by dp/dT < at intermediate T 
and low n, arising from the same temperature depen- 
dent disorder effect which produces the effective metal- 
lic phase [dp/dT > 0) at low T and high n. We find 
that for a range of densities below n c our calculated re- 
sistivity, p(T,n), as a function of temperature (T) and 
density (n), exhibits striking qualitative similarity to the 
experimentally observed p{T,n < n c ) [6-10], where the 
resistivity increases with decreasing T (i.e. dp/dT < 0) 
in a rather slowly rising approximate power law manner 
in contrast to the thermal exponential (either Mott vari- 
able range hopping or activated) behavior expected of 
a true T = insulator. The strongly insulating phase 
(which is beyond the scope of our theory) may occur at 
still lower densities where p(T) would rise exponentially 
with decreasing temperature. Our theory applies to all 
2D systems where the 2D MIT has been reported, and 
most significantly, in some systems (such as high mobil- 
ity 2D n- and p-GaAs) the observed 2D insulating phase 
seems to be entirely the novel low density effective in- 
sulating phase identified in this paper as the exponen- 
tial insulating behavior seems not to manifest itself in 
high-quality GaAs systems except perhaps at the lowest 
densities and temperatures. An important salient fea- 
ture of this low density novel phase (with dp/dT < 0) is 
that it is purely a temperature-induced crossover behav- 



1 



ior, and as such the critical density n c at which the 2D 
system crosses over to this phase (as n is decreased from 
above n c to below) is completely non-universal, being a 
function of temperature itself, i.e. n c = n c (T), with the 
effective critical density n c progressively decreasing with 
decreasing carrier temperature. The experimentally de- 
termined "critical" density is then n c (T = T m i n ) where 
T m in(~ 50 — 300mK depending on the experiment) is 
the effective minimum temperature to which the 2D car- 
rier system (and not the background bath) can be cooled 
down (typically T min 3> Tbath in 2D semiconductor sys- 
tems for low values of Tbath)- Such a non-universal be- 
havior for the critical density defining 2D MIT, namely 
that n c = n c (T) with the observed effective n c decreasing 
with decreasing temperature, has been pointed out ex- 
perimentally [9,10]. In particular, our theory specifically 
rules out a sharp (temperature independent) separatrix 
delineating the effective metallic (dp/dT > 0) and the 
effective insulating (dp/dT < 0) phase in the 2D MIT 
phenomenon. We emphasize that the experimental evi- 
dence for a sharp separatrix in the 2D MIT is extremely 
sparse, and most 2D MIT experiments observe a tem- 
perature dependent crossover density n c (T) separating 
dp/dT > for n > n c and dp/dT < for n < n c in 
agreement with our theory. 

We refer to this novel (low n/intermediate T) crossover 
"insulating" (dp/dT < 0) phase as an "apparent" insulat- 
ing phase in order to emphasize the fact that this phase 
is not a true insulating phase at T = 0, but is rather 
a low density (and "high" temperature) metallic phase 
where p(T) decreases with increasing temperature in a 
power law fashion. This newly identified effective insu- 
lating phase (i.e. dp/dT < 0) is only an apparent phase 
since lowering temperature further (at a fixed n < n c ) 
will eventually cause a temperature-induced re-entrance 
into the effective metallic phase at sufficiently low tem- 
peratures which, however, may not be accessible experi- 
mentally making the apparent insulating phase to appear 
to be a true insulating phase due to this low-T cut off. 
We emphasize that the hallmark of our proposed appar- 
ent insulating phase is the approximate power law tem- 
perature dependence of the resistivity in contrast to a 
true insulating phase where p(T) diverges exponentially 
with decreasing temperature. All localization effects, ei- 
ther weak or strong, are neglected in our theory. This is 
entirely consistent with our motivation of trying to un- 
derstand carrier transport in a "high temperature phase" 
where quantum interference effects should be small. The 
Boltzmann theory, described below, should be well-valid 
for this effective high-temperature semiclassical phase in- 
dependent of the actual value of the resistivity. 

We calculate the resistivity of a 2D carrier sys- 
tem scattering off random charged impurity centers 
with a screened Coulomb impurity scattering poten- 
tial. Within the Drude-Boltzmann semiclassical trans- 
port theory with the carrier-impurity interaction treated 
in the standard ensemble averaged Born approximation, 



we can express the 2D resistivity p as [11] 
p^ 1 = ne 2 (T)/m 7 



(1) 



where m is the carrier effective mass assuming a parabolic 
band effective mass approximation for the semiconductor 
band structure with n being the 2D carrier density. The 
thermally averaged scattering time (t) can be written, 
within the Drude-Boltzmann semiclassical theory (and 
the standard relaxation time approximation) , as 



(r) = (t)e = 



JdEr(E)E(-§l) 



(2) 



where (• • -)e indicates an energy averaging over the car- 
rier thermal distribution function f(E) = [l + exp{P(E — 
m)}] 1 with p = p(T,n) as the finite temperature chem- 
ical potential ((3 = l/fegT) given by p — ^ln[— 1 + 
exp((3Ep)] for a 2D system (note that p — Ep at T = 0, 
but p at finite T could be substantially different and is 
negative for T > 1.47V). The central quantity of interest 
is therefore the energy dependent transport relaxation 
time (i.e. the scattering time) t(E), which in the Born 
approximation is given by 
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where e k = H 2 k 2 /2m is the effective mass band energy of 
the 2D carriers with k = |k| as the 2D wave vector, and 
Ni is the three dimensional volume density of the random 
charged impurity centers of the a-th kind (we allow, in 
general, of different kinds of quenched disorder in the 
problem arising from the distributions of different types 
of random charged impurity centers in the 2D system) 
- we take z direction to be the confinement direction 
perpendicular to the relevant 2D plane of confinement of 
the carrier system. In Eq. (3) Uq*\z) is the 2D (in the 
x-y plane) Fourier transform of the 3D screened electron- 
impurity interaction in the system defined to be 
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where the effective dielectric function e(q) is given by 



2ttp 2 

e(q) = 1 - -^/(g)n(g), 



(5) 



where H(q) is the finite temperature static 2D noninter- 
acting polarizability function and the background 
static lattice dielectric constant. (We have carried out 
calculations also by including several different local field 
corrections in the 2D polarizability, but our qualita- 
tive results remain unaffected and therefore we will only 
present here RPA screening results where n is the 2D 
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non-interacting static polarizability.) In Eqs. (4) and 
(5) F(q, Zi) and f(q) are respectively the 2D subband 
form-factors for electron-impurity and electron-electron 
interactions, associated with the finite width of carrier 
wavefunctions in the z direction. These form factors 
are essential for quantitative accuracy, most particularly 
at the very low carrier densities of our interest where 
the quasi-2D electron layer spreads considerably in the 
transverse direction making the strict 2D approximation 
highly unreliable. We calculate the form factors by using 
the Fang-Howard-Stern variational wavefunctions for the 
quasi-2D layer [11]. For more details on the theory see 
refs. [11] and [12]. 

We have carried out extensive calculations of dc trans- 
port, based on the above-described Boltzmann theory, in 
several different 2D systems of current interest (i.e. sys- 
tems where an apparent MIT, as reflected in the change 
of sign in dp/dT at some low density, has been experimen- 
tally observed) at low temperatures (T = 50mK — 5K) 
and densities (n around, n ~ n c , the reported critical 
densities being n c w 10 11 , 10 10 , 10 9 cto~ 2 for typical 
high quality n-Si-MOS, p-GaAs, and n-GaAs systems re- 
spectively). Before presenting our numerical results for 
p{T,n) in Figs. 1 - 3, we first discuss qualitatively the 
theoretical results expected on the basis of Eqs. (1) - (5) 
in various temperature and density regimes. The dimen- 
sionless temperature, T/Tp ~ n _1 , plays an important 
role in our calculated temperature and density depen- 
dence 2D resistivity p(T, n). 

First, a straightforward asymptotic expansion of the 
Boltzmann transport equation shows that, for T/Tp <C 1, 
the resistivity can be written as p = po + Ap(T), where 
Po = p(T = 0) and Ap is given at very low temperature 
by [12,13] 

^-M£K(£f -*(£)"■ 

(6) 

where A\(> 0), Ay,/ 2 (> 0) arise from the thermal smear- 
ing of the electronic polarizability at 2kp wave vector, 
and B 2 (> 0) arises from the thermal energy averaging 
over the Fermi surface. Note that the first two screening 
terms in Eq. (6) are both positive, indicating an increas- 
ing resistivity with increasing temperature (the so-called 
2D metallic behavior) due to the suppression of screen- 
ing at finite temperatures whereas the third term, arising 
from the thermal energy averaging in Eq. (2), is a nega- 
tive correction to resistivity (i.e. a positive correction to 
conductivity) due to the increased thermal velocity of the 
carriers at finite temperatures. The expansions for A\, 
A 3 / 2 7 and B 2 have earlier been derived in the literature 
[5,12,13]: 

A 3/2 = 2.646 ( 1 + , (8) 



B 2 = —p(p+l), (9) 

where q$ = qTF/2kp, f = f(q = 2kp) is the 2D subband 
form-factor (defining the deviation from the strict two 
dimensionality in the layer) at wave vector 2/cf, and p 
(<~ 1) gives the energy dependence of the scattering time 
in Eq. (3), i.e. t(E) ~ E p for E <~ Ep. Three impor- 
tant feature of the asymptotic expansion given in Eq. (6) 
need to be emphasized: (1) The A and the B terms in Eq. 
(6) arise from totally different physics, and therefore we 
have neglected a possible A 2 (T/Tp) 2 term in the screen- 
ing expansion in Eq. (6) which may very well modify the 
B 2 term to B 2 = (B 2 — A 2 ). We use Eq. (6) only for 
the purpose of qualitative discussion, and all our results 
are based on the full numerical calculation of the semi- 
classical transport equations (Eqs. (1) - (5)) making any 
qualitative difference between B 2 and B' 2 irrelevant for 
our results. (2) There are higher-order interaction cor- 
rections for all the terms in Eq. (6), but these corrections 
are quantitatively unknown and are therefore left out of 
our discussion completely. (3) The asymptotic expansion 
only holds at very low temperature (T/Tp -c 1) which 
may not be of relevance to the 2D MIT phenomena [12]. 
Details on the low temperature transport behavior can 
be found in our recent work [12]. Here we are interested 
in the "high temperature" resistivity, which is manifestly 
not given by Eq. (6). 

It is possible to obtain the asymptotic high temper- 
ature result for the classical limit (T 3> Tp) of the 2D 
resistivity defined through Eqs. (1) - (6), and the result 
is [4] 

p(T»7»~(7>/T), (10) 

indicating that the "high temperature" resistivity de- 
creases with increasing temperature in an apparently 
"insulating-like" manner (i.e. dp/dT < for T ^> Tp). 
This can be rather easily understood physically based on 
the fact that at high temperatures the main effect arises 
from the thermal increase in the average carrier velocity 
which obviously leads to a decreasing resistivity. Thus, 
the linear power law resistivity, p ~ T^ 1 , in Eq. (10) 
arises from the energy averaging in Eq. (2). 

Even before carrying out the full numerical calcula- 
tions it is obvious that any smooth interpolation between 
the "low temperature" behavior of Eq. (6) and the "high 
temperature" behavior of Eq. (10) would exhibit non- 
monotonic temperature dependent resistivity [4] associ- 
ated with this quantum-classical crossover, which has 
been seen [6-10] experimentally. This leads to the fol- 
lowing generic temperature dependence of p(T/Tp): At 
very low temperatures p is linear in T/Tp with dp/dT 
(> 0) being determined by the absolute value of the co- 
efficient A x ; with increasing T/Tp (but with T/T F < 1) 
it is possible (but not necessarily essential) for dp/dT to 
change sign becoming negative, dp/dT < 0, if the B 2 
term can overcome the effects of the first two positive 
terms within the constraint T/Tp < 1. This obviously 
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depends on the relative values of coefficients A\, A 3 / 2 , 
and £?2- In Si (100) inversion layer and in p-GaAs sys- 
tem screening effects are strong leading to A\, A 3 / 2 3> B2 
so that at low T/Tp the metallic resistivity, dp/dT > 0, 
dominates. In n-GaAs 2D system, on the other hand, 
screening effects are weak (since the density of states 
is small) and for modest values of T/Tf, the Bi term 
is larger than the first two terms, particularly at high 
densities where 2kp ^> qTF, leading to dp/dT < for 
moderate values of T/Tp. (In fact, this is also expected 
in highly disordered systems where screening effects may 
be cut off by impurity scattering at a characteristic tem- 
perature To which may not be too low.) The message 
from the analytical results of Eqs. (6) - (9) is therefore 
clear: A 2D system would always exhibit "metallic" re- 
sistivity (defined as dp/dT > 0) for the lowest value of 
the dimensionless temperature T/Tp (ignoring all weak 
localization effects which we neglect throughout this pa- 
per) crossing over to an apparent "insulating" (defined as 
dp/dT < 0, but not with an exponential temperature de- 
pendence) resistivity at some intermediate temperature, 
Across j where "metallic" (dp/dT > 0) screening effects are 
overwhelmed by the "insulating" (dp/dT < 0) thermal 
averaging effects. If the carrier density of the 2D system 
is effectively low so that T cross is lower than the lowest 
achievable carrier temperature (T m ; n ) in the system, then 
the 2D system will manifest an apparent insulating phase 
(dp/dT < 0) down to the lowest experimental tempera- 
ture T (> T cross ) since the low temperatures (T < T cross ) 
where the screening effect dominates for a crossover to 
the metallic behavior (dp/dT > 0) cannot be achieved. 
We note that it is possible in many situations, depending 
on carrier densities and system parameters (e.g. effective 
mass, type of impurity scattering, etc.), for T cross < Tp. 
Remembering that Tp w 0.7n K (Si and p-GaAs) and 
4n K (n-GaAs) where n is the 2D carrier density mea- 
sured in the units of 10 10 cm -2 , we conclude that Tcross 
could be as low as lOOmK (in some situations) for 2D 
carrier densities around 10 10 cto~ 2 for n-Si (and p-GaAs) 
and 10 9 cto~ 2 for n-GaAs systems. A characteristic fun- 
damental feature for the crossover we are discussing here 
is the strong dependence of the crossover temperature on 
the carrier density with T cross (n) increasing strongly with 
carrier density. The other important feature is the non- 
existence of the apparent insulating (dp/dT < 0) phase 
at low enough temperature (T < T cvoss ), i.e. if the carrier 
temperature of a 2D system can be decreased arbitrar- 
ily then (leaving out all localization effects) the system 
would always become metallic (dp/dT > 0) at low enough 
temperatures (defined as T <C T cross ). 

Our numerical transport results shown in Figs. 1-3 
for p-GaAs and n-GaAs respectively bear out the theo- 
retical considerations outlined above. In particular, the 
calculated resistivity is always metallic at the lowest tem- 
peratures, manifesting a clear crossover to an insulating 
behavior at T CTOSS (n) with a change of sign in dp/dT (i.e. 
dp/dT > for T < T CIOSS and dp/dT < for T > r croS s) 
as described above. We note that T cross could be quite 
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FIG. 1. (a) Calculated 2D resistivities for p-GaAs sys- 
tem for various hole densities p =0.1, 0.2, 0.5, 1.0, 2.0, 
3.0xl0 10 cm~ 2 (from the top) as a function of temperature. 
In this figure we use the charged bulk impurities inside the 
quasi-2D systems, (b) Calculated 2D resistivities for n-GaAs 
system for various electron densities n =0.4, 0.5, 1.0, 2.0, 3.0, 
4.0, 5. Ox 10 9 cm -2 (from the top) with the charged interface 
impurities as a function of temperature. (T D = is used for 
these results.) 

low for low values of carrier density, and therefore it 
is possible, in principle, for the apparent insulating be- 
havior to arise entirely from this crossover effect. In 
particular, in Fig. 1 we choose (rather arbitrarily) 
the lowest possible achievable carrier temperature to be 
300mK — 500mK (for the purpose of illustration) and 
show what a typical 2D MIT transport data set (in the 
T = 300mK— 3K range) would look like based entirely on 
this finite temperature crossover phenomenon. The over- 
all qualitative consistency between our results and many 
of the existing 2D MIT experimental results in the litera- 
ture suggests that the simple (and physically appealing) 
crossover scenario presented in this paper cannot be en- 
tirely ruled out as the cause for the 2D MIT behavior 
observed experimentally at least in some situations. We 
note a quantitative discrepancy between our theory and 
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experiment: in general, we find X^ross in our theory to 
be higher than the corresponding crossover temperature 
in the experiment. (This discrepancy can be corrected 
by introducing a cut-off To in the screening as discussed 
later in this paper.) 

It is essential to emphasize two points: (1) The effec- 
tive insulating phase is only an apparent insulating phase 
by virtue of the insulator being defined through the nega- 
tive temperature derivative of resistivity, i.e. dp/ dT < 0. 
This apparent insulating phase is really a metallic (i.e. 
extended or delocalized electron wavefunction) electron 
liquid phase (interacting with random charged impurity 
centers) exhibiting dp/ dT < down to a crossover tem- 
perature scale (T > T cross ) with the lowest tempera- 
ture (T <C r cr oss) phase being always a true metal with 
dp/dT > 0. (2) The temperature dependence of p(T) 
in this apparently insulating phase simulates a power 
law (with p ~ 1/T in the T > Tf high-temperature 
region) rather than the exponential temperature depen- 
dence typical of a true insulating phase. The finite tem- 
perature crossover "insulating" phase (i.e. dp/dT < 0) is, 
by definition, a "high" temperature phase in our Drudc- 
Boltzmann transport model with the true "low" temper- 
ature (T — ► 0) phase in our theory (which manifestly 
excludes weak and strong localization effects) being al- 
ways a metallic (i.e. dp/dT > 0) phase. The crucial 
point of physics is, however, the fact that experimentally 
there is always a low temperature cut-off, T m ; n , below 
which carriers simply cannot be cooled down, and for 
Tcross < T m j n , the 2D apparent "insulating" phase dis- 
cussed in this paper would, for all particular purposes, 
behave like a 'real' insulating phase. The "apparent" 2D 
metal-insulator transition in this crossover scenario will 
therefore appear to occur at the density n c defined by 

Tcross(^c) ~ ^min- 

In comparing our theoretical results with the existing 
2D MIT data in the experimental literature, we men- 
tion that a large number of reported data look quali- 
tatively similar to the results shown in Figs. 1 - 3 of 
this paper. In particular, the following features of our 
theoretical results are often observed experimentally: (1) 
the low-density-low temperature 2D insulating phase (for 
n < n c ) is often seen to have an approximate power law 
temperature dependence resistivity; (2) the low tempera- 
ture extrapolated (T — ► 0) resistivity does not seem to be 
divergent in many cases; (3) the experimentally observed 
critical density n c has been reported to be a function of 
the lowest measurement temperature with n c (T) decreas- 
ing monotonically with decreasing temperature. For low 
enough carrier densities, most 2D systems should even- 
tually enter a strongly localized transport regime where 
the low-temperature resistivity diverges in an exponential 
manner, but in high quality (i.e. low disorder) systems 
of interest in the 2D MIT phenomenon such exponen- 
tially divergent resistivity seems to be always preceded 
by a density regime (below n c ) where p(T) increases 
smoothly in a power law manner with decreasing tem- 
perature. This intermediate density regime of apparent 



insulating behavior is most strongly manifested in 2D p- 
GaAs and n-GaAs systems, but it has also been seen in 
n-Si MOS systems. A possible reason for the prevalence 
of the apparent insulating phase in 2D GaAs systems 
could be the much higher quality (i.e. much lower dis- 
order) for GaAs-based 2D systems compared with 2D Si 
structure, making it possible for the system not to mani- 
fest localization down to rather low carrier densities, thus 
enabling the crossover insulating phase to show up more 
easily without being preempted by localization. 

In discussing experimental data, it is important to real- 
ize that the measured 2D resistivity p(T) typically shows 
a saturation at low enough temperature, i.e. below some 
saturation temperature (T s ~ 50 — 300mK, depending 
on the experiment) p(T < T s ) becomes essentially a con- 
stant completely independent of temperature. Although 
an underlying (unknown) fundamental cause for this re- 
sistivity saturation behavior cannot be definitively ruled 
out, the saturation is generally believed to arise from 
a saturation in the carrier temperature due to the in- 
evitable carrier heating problems in semiconductors at 
low temperatures and densities. This possibility of car- 
rier temperature saturation at low temperatures (where 
further cryogenic cooling only lowers the temperature of 
the surrounding bath, not the 2D electrons themselves) 
is further supported by the fact that the resistivity satu- 
ration is also seen in the insulating system at low enough 
carrier temperatures. We have completely ignored the 
low temperature resistivity saturation problem in our 
analysis, assuming, perhaps somewhat uncritically, but 
in agreement with all existing theoretical analyses in the 
subject, that the resistivity saturation arises from tem- 
perature saturation, i.e. T s = T min . Within the screening 
theory, in fact, thermal suppression of screening will be 
cut off at temperatures (T < To) lower than the temper- 
ature scale of the disorder scattering (with To > h/ksT), 
and we have taken this disorder broadening effect approx- 
imately into account in the results presented in Figs. 2 
and 3 of this paper. Another comment in the context of 
2D MIT experiments is the issue of the so-called separa- 
trix, where the critical density n c separating the metallic 
and the insulating phase is claimed to manifest a com- 
pletely temperature independent resistivity sharply de- 
lineating the metallic phase (dp/dT > 0) from the in- 
sulating phase (dp/dT < 0) and thus defining a quan- 
tum phase transition. Such a sharp separatrix cannot 
be explained at all within our Drudc-Boltzmann theory 
since all electronic states are by definition extended in 
our semiclassical theory. As has been emphasized above, 
the apparent critical density n c in our theory is strongly 
temperature dependent with n c delineating dp/dT > 
(for n > n c ) and dp/dT < (for n < n c ) decreasing with 
the lowest measurement temperature T m ; n with n c — > 
as T m i n — > since a true insulator does not exist in our 
theory. An experimental separatrix defining the 2D MIT 
has, however, been rarely reported, and in that sense 
our theory is consistent with a large body of the existing 
2D MIT experimental data. We do emphasize, however, 
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that at low enough densities all real 2D systems should 
eventually exhibit exponential localization behavior (at 
densities below the "power law" apparent localization be- 
ing discussed in this paper) which is beyond the scope of 
our work. We believe that Si MOS systems, typically 
being more disordered than 2D GaAs systems, exhibit 
much more of the exponential-type true insulating be- 
havior than the power-law apparent insulating behavior 
being discussed in this paper. 

In a very recent paper Noh et al. [7] have reported an 
anomalous power law insulating behavior in a 2D p-GaAs 
systems at low densities (and intermediate temperature), 
which we believe to be the apparent insulating phase in- 
duced by quantum-classical crossover being predicted in 
this paper. The r s values in the samples of rcf. [7] are 
extremely high (r s ~ 44 — 80), bringing up the possibil- 
ity of Wigner crystal (perhaps in a strongly correlated 
molten phase [14]) physics playing a role in this exper- 
iment. We believe that it is unlikely that Wigner crys- 
tal (WC) physics has much to do with this anomalous 
insulating behavior. First, the experimental tempera- 
ture and density regime where the anomalous insulat- 
ing behavior is observed is well outside the WC phase 
boundary (see our Fig. 4 in this paper). Second, the 
transport data plotted as a(T) = p^ 1 for various densi- 
ties in ref. [7] show absolutely no indication of a critical 
density separating metallic and insulating phases — the 
data at all densities look essentially the same at "high 
temperatures" where a(T) <~ T whereas at "low tem- 
peratures" the relatively high density conductivity shows 
an upward bending below some density-dependent tem- 
perature T cross (n) (i.e. dp/dT changes sign from being 
insulating-like, dp/dT < 0, at higher temperatures to be- 
ing metallic-like, dp/dT > 0, at lower temperatures). We 
emphasize that for T > T cross (n) there is absolutely no 
difference in the observed qualitative behavior of cr(T) at 
different densities. Since T cross (n) decreases with carrier 
density n, it is impossible to rule out the possibility that 
the experimental low density a(T) plots would actually 
also bend upward at temperatures lower than the lowest 
hole temperature (~ 60mK) achieved in these experi- 
ments. In fact the experimental data of Noh et al. arc 
entirely consistent with this scenario. For example, if the 
experimental temperature cut-off is taken to be 200mK 
in the data of Noh et al, then all the plots shown in 
ref. [7] would appear to be "insulating" (upto the hole 
density of 5 x 10 9 cm~ 2 ) since all the experimental data 
exhibit dp/dT < (i.e. da/dT > 0) down to 200mK. 

To emphasize the similarity between the experimental 
data of Noh et al. and our predicted crossover appar- 
ent insulating phase we show in Fig. 2 our calculated 
low density transport results for the Noh et al. sam- 
ple. Remembering that phonon scattering (neglected in 
our theory) becomes important for GaAs holes already 
at 500 mK [6], the results of Fig. 2 are remarkably simi- 
lar to the experimental data presented in rcf. [7] . This is 
particularly true for Fig. 2(b), where we include a Dingle 
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FIG. 2. Calculated conductivities for p-GaAs system [7] 
for various hole densities p =5.0, 4.0, 3.0, 2.3, 1.5xl0 9 cm -2 
(from the top) with the remote charged impurities. In (a) we 
calculate conductivities without level broadening in the polar- 
izability function, and in (b) we include the level broadening 
in the screening, following ref. [14], with Td =0.2, 0.25, 0.3, 
0.4, 0.5K (from the top). The vertical bars indicate the Fermi 
temperature T F . 

temperature cut-off in the screening function in order 
to approximately simulate the effect of collisional broad- 
ening (due to impurity scattering). The basic idea [15] 
is to include in the finite wave vector static polarizabil- 
ity the level broadening (parametrized by a Dingle tem- 
perature Td = Y/kg where V is the collisional damp- 
ing) arising from the impurity scattering. As discussed 
earlier in the literature [15] such a collisional damping 
term, parametrized by the level broadening parameter 
Td, acts to suppress the thermal effects on screening at 
low temperatures, T <C Td, while leaving the screen- 
ing function essentially unaffected at high temperatures 
T > Td- Since the level broadening parameter in the 
screening, Td, is unknown (within our approximation 
scheme), we do not attach any great quantitative sig- 
nificance to our approximation except to note that this 
zcroth order scheme of incorporating collisional damping 
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in screening is physically meaningful since impurity scat- 
tering should tend to suppress finite wave vector static 
screening of an electron gas. We emphasize that the pa- 
rameter Td used in our theory should not be construed 
to be Td ~ ft/fc^r (which would be far too small), but 
should be considered an effective level broadening param- 
eter which cuts off the strong temperature dependence of 
screening for T <Td, with Td being a parameter of the 
theory (increasing with decreasing conductivity) . 

It is therefore worth while to emphasize that our results 
(with Td ^ in screening) shown in Fig. 2(b) are quali- 
tatively (and even semiquantitatively) in good agreement 
with the approximately linearly rising, a oc T, conduc- 
tivity at higher temperatures found by Noh et al. In fact, 
even our Td = static RPA screening transport results 
are in qualitative agreement with experiment except for 
the fact that the crossover temperature T cross for da/dT 
to change sign is consistently much higher in Fig. 2(a) 
with Td =0 than in our Fig. 2(b) with Td ^ or in the 
experiment. The basic behavior of an approximately lin- 
ear temperature dependence of conductivity, a cx T, for 
T > T CIOSS (n) with T cross (n) decreasing with decreasing 
density (and da/dT < 0, i.e. er(T) increasing with de- 
creasing temperature for T < T cross (n)) applies to Figs. 
2(a) and (b) and to the experimental data. Inclusion 
of phonon scattering (neglected in our calculations here) 
is likely to further improve the good agreement between 
experiment and theory. 

In the context of comparison with the experimental 
data we show in Fig. 3 our calculated low temperature 
transport results for the 2D GaAs electron system cor- 
responding to the recent measurement of Lilly et al. [8] 
carried out in a high quality low density n-GaAs gated 
heterostructure system. We have kept the level broaden- 
ing parameter Td — 0.75K fixed in the results shown in 
Fig. 3 and have shown the behavior of both the resistiv- 
ity p(T) and conductivity a(T) for several low densities. 
Again, phonon scattering effects, possibly of some impor- 
tance here for T > IK [8], have been neglected. These 
theoretical results are again in striking qualitative and 
semiquantitative agreement with the experimental data 
of Lilly et al. [8] demonstrating that the experimentally 
observed anomalous insulating phase in high-quality (i.e. 
low-disorder) 2D semiconductor systems may very well 
be the apparent crossover insulating phase (which would 
manifest metallic behavior for T < T cross (n) except that 
T cross may be too low to be experimentally accessible at 
low carrier densities) predicted by our theory. 

The striking qualitative agreement between our the- 
ory and experiment for the anomalous insulating phase 
encourages us to predict the following scenario for high 
quality GaAs-based 2D (both electron and hole) sys- 
tems. As 2D carrier density is lowered the high tem- 
perature crossover insulating phase (T > T cross ) extends 
to lower temperatures as T CIOSS (n) goes down, and even- 
tually when Tcross goes below the lowest accessible exper- 
imental temperature, the system appears to be insulating 
with an anomalous (i.e. power law) temperature 
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FIG. 3. Calculated (a) resistivities and (b) conductivities 
for n-GaAs system for various electron densities n =1.6, 2.0, 
2.3, 2. 9xl0 9 cm -2 with the remote charged impurities and a 
fixed Dingle temperature To = 0.75K as a function of tem- 
perature. The vertical bars indicate the Fermi temperature 
T F . 

dependence characterizing the apparent insulating phase. 
Since our predicted scenario is in excellent qualitative 
agreement with the available experimental data in GaAs- 
based high-mobility 2D holes and electrons, we propose 
the following direct experimental test for its verification 
(or falsification): One should lower the temperature to 
check whether the observed anomalous temperature de- 
pendence of the putative insulating phase remains insu- 
lating or bends upward (i.e. da/dT goes through zero 
changing its sign from positive to negative as it does 
at higher densities) at lower temperatures. Since low- 
ering carrier temperature arbitrarily is typically difficult 
to achieve in 2D semiconductor systems, a relatively easy 
way of testing our proposed scenario will be to directly 
plot the experimentally measured da/dT (or dp/dT) as a 
function of temperature at low densities in the anomalous 
insulating phase. If the magnitude of this derivative is 
decreasing with decreasing temperature indicating that 
a (or p) may be approaching a minimum (or maximum) , 
then this is a rather strong verification of our proposed 
scenario. We mention in this context that the experi- 
mental data of ref. [7] and [8] are qualitatively consistent 
with our proposed scenario of a decreasing magnitude of 
da/dT as T decreases. We believe that there are simply 
two possibilities (which are consistent with our predic- 
tion of the anomalous power law insulating phase being 



7 



a crossover phase): Either the magnitude of dp/dT will 
decrease indicating a transition to an effective metallic 
phase or dp/dT will diverge exponentially indicating a 
true insulating phase. We suggest that experiments be 
carried out in the anomalous power law insulating phase 
(with careful measurements of dp/dT as a function of T 
and n) to verify our prediction. 

Finally, we discuss the possibility of WC physics play- 
ing a role in the observed anomalous insulating behavior 
of refs. [7,8]. In Fig. 4 we show our calculated approxi- 
mate WC (T — n) phase diagram for both 2D electrons 
and holes in GaAs, combining both the T = quan- 
tum Wigner crystallization at low carrier densities and 
the classical Wigner crystallization. The calculation of 
this approximate WC phase diagram follows our recent 
work in ref. [16] where we show how one can implement 
an interpolation scheme to obtain the full phase diagram 
in the density-temperature plot by combining the known 
classical [17] and quantum [18] WC limits. We also show 
as shaded regimes the experimental samples correspond- 
ing to refs. [7] and [8] in Fig. 3(a) and (b) respectively. 
Both the shaded regions are well outside the WC phase in 
the phase diagram indicating that WC physics is unlikely 
to be playing a dominant role here. The more important 
reason for ruling out the WC scenario is, however, ex- 
perimental — the experimental data of Noh at al. [7] is 
completely smooth and the conductivity behavior in the 
"metallic" (higher density) and the "insulating" (lower 
density) phase is identical at higher temperatures, the 
only difference being the upward curvature in er(T) in 
the higher density data which moves to lower temper- 
atures in the lower density plots and may have simply 
moved to temperatures below the measurement temper- 
atures at the lowest densities making it experimentally 
invisible. 

In summary, we have identified an apparent "high tem- 
perature" insulating phase in 2D systems at intermediate 
carrier densities where the conductivity increases with 
temperature in a power law fashion (approximately lin- 
early). This apparent or effective insulating phase is a 
temperature-induced quantum-classical crossover phase 
which exists only for T > T cross (n) where T cross decreases 
with decreasing carrier density. We have shown that, al- 
though typically T cross ~ Tp (oc n), it is possible for 
Across to be substantially below the Fermi temperature 
if collisional broadening effects are included in the car- 
rier screening function. We have shown that our calcu- 
lated transport behavior in this effective insulating phase 
agrees qualitatively very well with recent observations in 
high-quality 2D electrons and holes in GaAs heterostruc- 
tures, We have not discussed the 2D Si-MOS system in 
this paper because the disorder effects in Si MOSFETs 
are substantially higher (mobilities are typically factors 
of 10 — 1000 lower in Si MOS systems than in GaAs sys- 
tems), and the Si MOS systems make transitions to the 
strongly localized phase (with exponential temperature 
dependence) much sooner — this intermediate crossover 
insulating phase therefore exists only in a very narrow 
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FIG. 4. Calculated 2D hole/electron liquid - Wigner crys- 
tal phase diagram in the density - temperature (n, T) plot for 
(a) p-GaAs and (b) for n-GaAs systems. The shaded regions 
indicate the experimental samples corresponding to ref. [7] for 
p-GaAs and ref. [8] for n-GaAs. Here Tp is the Fermi tem- 
perature and T = (V) I '(T) is the ratio of the classical mean 
potential energy to the mean classical kinetic energy. The 
classical WC - electron liquid (first order) phase transition line 
is defined by T — T c , which is found from numerical (molec- 
ular dynamics) simulations [17] to be F c « 120. At T — 
a 2D system is expected to become a quantum WC phase at 
r s > 37 [18]. The heavy solid line, indicating the liquid-solid 
phase boundary, is obtained from the interpolating scheme 
described in ref. [16]. The region above (below) the T = TV 
line is the classical (quantum) region with CWC (QWC) de- 
noting the classical (quantum) Wigner crystal phase. 

range of densities and temperatures in Si MOS systems. 

The weakest empirical point in our quantitative com- 
parison with the experimental data is the fact that 
our theoretically obtained T^ oss (n) seems to be consis- 
tently higher than the experimental crossover tempera- 
ture T c e r o SS (n) where T cross at a particular density could 



for example be defined by the condition da/dT = 
at T — T cross . We have tentatively fixed this problem 
in this work by introducing the collisional damping in- 
duced screening cut-off through the Dingle temperature 
parameter Tp,. The other possibility is to take the real 
Fermi temperature Tp to be smaller than the nominal 
Tp (~ n/m) given by the carrier density n and the band 
effective mass m which could happen if the actual effec- 
tive mass is larger than the band mass and/or the actual 
free carrier density is smaller (e.g. due to trapping of car- 
riers). While in the GaAs-based high mobility systems 
studied in this paper, the discrepancy between T*^ oss (n) 
and T^oss ( n ) 1S relatively small even for the To = the- 
ory, this discrepancy is very large in Si MOS based 2D 
systems where disorder is high. A straightforward appli- 
cation of our Tp, = theory for Si MOS systems indicates 
T cr oss > IK for the n = 5x 10 10 -10 n cm~ 2 density range 
of relevance to Si 2D systems. To reduce T croS s substan- 
tially one needs to introduce large values of T D where 
the theory is not meaningful. The other possibility is 
that the effective mass in Si may be enhanced from its 
band value and/or the effective carrier density may be 
substantially lower than the nominal n. It is therefore 
interesting to note that there is no clear cut signature of 
this anomalous insulating phase in 2D Si MOS systems 
where the insulating phase, for T < IK or so, indicates 
exponentially rising resistivity with decreasing tempera- 
ture. This is consistent with our finding of rather large 
T^ oss in Si MOS systems implying that this apparent 
insulating phase in Si systems remains only a high tem- 
perature phase since the electrons can always be cooled 
below T cross in the 2D Si samples. 

We note finally that the validity of the semiclassi- 
cal Drudc-Boltzmann theory for studying the anomalous 
insulating phase is not a crucial issue in this context. 
The reason is that Boltzmann theory, being semiclassi- 
cal, becomes increasingly valid at higher temperatures, 
and therefore, should be a reasonable qualitative descrip- 
tion for the "high-temperature" effective insulating phase 
being discussed here. The precise value of the resistiv- 
ity p and the corresponding localization parameter fcpZ 
(where I is the transport mean- free path), i.e. whether p 
is smaller than the quantum of resistance h/e 2 or kpl is 
larger than unity, are important issues for T = consid- 
erations (where quantum localization effects are impor- 
tant) but not for finite T/Tp (« 0.1 — I) case which is 
of interest here. In fact, within our scmiclassical Boltz- 
mann transport theory the behavior shown in Figs. I - 
3 are generic behavior determined entirely by n and T, 
completely independent of the actual resistivity values of 
the 2D system, i.e. one should think of the ordinate (p 
or a in our figures) to have completely arbitrary units 
— the actual values of p (= er -1 ) in our theoretical re- 
sults being determined by the unknown charged impurity 
density (JVj) in the system which we use as an adjustable 
parameter to set the resistivity scale. It is more appro- 
priate to think of ordinate at Figs. 1 and 2 as p/po or 
ct/cto where p n or <r are the T = Drude values of resis- 



tivity or conductivity. These considerations make it clear 
that our qualitative results are valid (at finite T/Tp) ir- 
respective of the actual sample resistivities as long as one 
can neglect quantum localization effects which seems to 
be the case for the GaAs experiments in refs. [7] and [8], 
but perhaps not for Si MOS samples. 

In our conclusion we discuss critically the approxi- 
mations of our theory which may limit its applicability 
to real 2D systems. Our most important approxima- 
tions are the Drude-Boltzmann scattering theory (due 
to charged impurity scattering) and RPA screening by 
the electron liquid. Both of these approximations are 
simplistic at the low carrier density phenomena of in- 
terest here. But there are good reasons to believe that 
our theory is sound in explaining the qualitative behavior 
of 2D semiconductor systems at densities and tempera- 
tures above the true localization regime. First, recent 
systematic diagrammatic calculations [5] of higher order 
interaction corrections to the 2D resistivity show that 
the basic picture of an effective 2D metallic behavior with 
dp/dT > (with a leading order linear temperature coef- 
ficient of resistivity) applies well in the ballistic transport 
regime (h/hsT <C T <C Tp) provided weak localization 
effects are negligible. Second, the results for p(T, n) ob- 
tained in the 2D metallic regime within the RPA-Drudc- 
Boltzmann theory agree well with existing experimen- 
tal results. Third, RPA becomes quantitatively accurate 
at high temperatures (and, our inclusion of local field 
correlation corrections going beyond RPA give qualita- 
tively very similar results). We therefore believe that our 
description of the high-temperature apparent insulating 
phase based on the semi-classical Drude-Boltzmann- RPA 
transport theory may have considerable theoretical valid- 
ity. 
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